Cluster magnetic fields from active galactic 

nuclei 

P.M. Sutter*, P.M. Ricker 4 -** and H.-Y. Yang f 

* Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801-3080 
^Department of Astronomy, University of Illinois at Urbana-Champaign, Urbana, IL 61801 
** National Center for Sup ere omputing Applications, University of Illinois at Urbana-Champaign, 

Urbana, IL 61801 

Abstract. Active galactic nuclei (AGN) found at the centers of clusters of galaxies are a possible 
source for weak cluster-wide magnetic fields. To evaluate this scenario, we present 3D adaptive 
mesh refinement MHD simulations of a cool-core cluster that include injection of kinetic, thermal, 
and magnetic energy via an AGN-powered jet. Using the MHD solver in FLASH 2, we compare 
several sub-resolution approaches that link the estimated accretion rate as measured on the simula- 
tion mesh to the accretion rate onto the central black hole and the resulting feedback. We examine 
the effects of magnetized outflows on the accretion history of the black hole and discuss the ability 
of these models to magnetize the cluster medium. 
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INTRODUCTION 

Faraday rotation measures of galaxy clusters suggest they contain magnetic fields with 
average strengths of ~ 1 jitG (e.g., yj]), but the origin of these fields is unknown. 
Possible seed fields, such as those generated by the Biermann battery mechanism, 
require amplification |0]. Accretion disks around active galactic nuclei (AGN) are a 
promising candidate amplification source, since AGN-launched jets require magnetic 
fields and the energetics of AGN-blown bubbles suggest that AGN can amplify the seed 
fields to the observed strengths. 

We would like to directly test this scenario in simulation, but since AGN accretion 
disks are typically smaller then ~ 100 AU, and cosmological simulations have resolu- 
tions on the order of ~ 1 kpc, we must include the amplification and jet launching as a 
subgrid model. We investigate two commonly-used AGN feedback models: one based 
on jets [3], the other on already-inflated bubbles [4!]. We link the feedback energies of 
these models to the strength of an injected magnetic field. We inject magnetic fields with 
the form described by [5], in which the field has both toroidal and poloidal components. 

ANALYTICAL AND NUMERICAL APPROACH 

We assume Bondi accretion: 



M Bo ndi = 4^G 2 m| H p/c s 3 , 



(1) 




where the sound speed c s and the density p are measured on the simulation mesh, and 
m BH is me black hole mass. Following [4], to compensate for under-resolving the actual 
accretion disk, we assume a constant multiple of the Bondi rate: 

M = aM B ondi (2) 

with a = 100. 

We follow [H] for modeling jet-based feedback where the energy injection rate is 

£ = £ F Mc 2 (l-l/M load )|n (3) 

Similarly, the momentum injection rate is P = ^/le^Mc^V , and the mass injection rate 
isM inj = M load M\ x ¥\. 

The window function x ¥, which provides a mapping onto the mesh, is 

n^ = ^T^P\-'-rj-)iT- (4) 

We cut off injection at z = /? e j and r = 2.6r e j. In the above, c is the speed of light, and we 
will define the injection region through r e j = 3.2 kpc and /? e j = 2.5 kpc. The injection 
region is oriented along the z-axis. We assume a jet mass loading factor of Mi oa d = 100 
and feedback efficiency of £p = 0.1. 

For bubble injection, as in [0], we have only thermal energy injection: 

E = £ F £ M Mc 2 : (5) 

where £p is the same as above and £m = 1.0. We distribute this energy uniformly in a 
sphere with radius determined by 

where we define the scalings by Rq = 43 kpc, Eq = 10 60 erg, p = lO 6 M kpc~ 3 , and 
dt is the timestep. These scalings ensure that a bubble in a typical cluster environment 
will have a realistic size. The bubbles are always centered on the black hole, and we 
only form bubbles when the black hole has increased its mass since the previous bubble 
formation by AMbh/^bh > 0.001. 

The injected magnetic field takes the form 



B r (r',z') = 2fl zVexp(-/ 2 -z /2 ) (7) 
B z (r',z') = 2£ (l-/ 2 )exp(-/ 2 -z /2 ) (8) 
= floa B /exp(-/ 2 -z /2 Y (9) 



where r 1 = a/x 2 +y 2 /ro and zf = z/tq. Here, r$ is l/2i? bub for bubbles and l/2i? e j for 
jets, and Ob is the ratio of of polodial to toroidal flux. We choose ocb = V^IO for an 
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FIGURE 1. Accretion rate as a fraction of the Eddington rate for magnetized and unmagnetized jets 
and bubbles. 

initially relaxed field, as suggested by [5] . We determine the scale Bq by giving half of 
the available feedback energy to the magnetic field. 

We performed three-dimensional simulations with an isolated cluster profile in a 2048 
kpc box using FLASH 2.5 fla], an adaptive mesh refinement (AMR) code. Both the jet 
and bubbles runs used a maximum resolution of 1 .0 kpc within the central 50 kpc region. 
The AGN began as a 1O 7 M black hole in the center of an NFW (170) gravitational 
potential. The cluster had a concentration of 6.5, scaling radius of 165 kpc, total mass 
of lO 14 h _1 M , and gas fraction of 0.12. We included cooling from [8] assuming 1/3 
solar metallicity. We allowed the cluster to relax for ~ 1 Gyr before activating cooling 
and feedback. 

RESULTS 

Figured] shows the accretion history of the black hole with both magnetized and unmag- 
netized jet- and bubble-based feedback until the cooling catastrophe occurs. We find that 
magnetizing bubbles does not greatly alter the accretion rate relative to unmagnetized 
injections, since bubbles occur infrequently. However, magnetizing jets greatly reduces 
the accretion rate. Here, the combination of an axial jet and a toroidal magnetic field 
prevents gas from accreting. 

We show the volume magnetized in Figure|2] Both jets and bubbles are able to weakly 
magnetize large volumes. These fields may be further amplified in realistic clusters by 
turbulence and merger shocks. However, significant (> 1 /iG) fields do not penetrate far 
from the cluster core. Finally, only the jet is able to produce greater than 20 /iG fields, 
but these quickly dissipate. The outflows begin to greatly enhance the magnetization of 
the cluster at the onset of the cooling catastrophe. 

With the parameter values used, we find that these feedback mechanisms do not 
provide enough heating to prevent the cooling catastrophe. We are investigating the new 
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FIGURE 2. Magnetized volume for various thresholds. The y-axis is the the cube root of the total 
volume that lies above each threshold. 

MHD solver in FLASH 3 ([@]) to determine whether or not this is a numerical effect. 

The evolution of our injected magnetic fields resembles that of |@], but they consider 
only pure magnetic injection from a single source with fixed energy input. Our energy 
input depends on the accretion rate, and we find much less available energy than the 
value they assume (~ 10 60 ergs). 
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